Method to characterize spinons as emergent elementary particles 
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We develop a technique to directly study spinons (emergent spin S — 1/2 particles) in quantum 
spin models in any number of dimensions. The size of a spinon wave packet and of a bound pair 
(a triplon) are defined in terms of wave-function overlaps that can be evaluated by quantum Monte 
Carlo simulations. We show that the same information is contained in the spin-spin correlation 
function as well. We illustrate the method in one dimension. We confirm that spinons are well defined 
particles (have exponentially localized wave packet) in a valence-bond-solid state, are marginally 
defined (with power-law shaped wave packet) in the standard Heisenberg critical state, and are not 
well defined in an ordered Neel state (achieved in one dimension using long-range interactions). 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.Mg, 75.10.Pq 



Spinons are emergent spin S — Xjl particles (fractional 
excitations) of quantum magnets [l|-|3| and potentially 
exist also in strongly-correlated electron systems such as 
the high-Tc cuprate superconductors [^j] . Their existence 
is well established in one-dimensional (ID) systems [liy , 
where they correspond to kinks and antikinks (solitons). 
In higher dimensions, gapped magnons ("triplons") can 
be viewed as bound states of spinons. Under some con- 
ditions, in spin liquid states |3| and at certain quantum- 
critical points [5], these spinons may become deconfined 
(unbound). Even in cases where the spinons are not 
completely deconfined, such as in a valence-bond-solid 
(VBS) state of a two-dimensional (2D) system close to a 
phase transition into the antiferromagnetic (Neel) state, 
the bound state can become very large [5|. The spinons 
can then be viewed as deconfined below the length-scale 
of the pair size, and above a corresponding (relatively 
low) energy scale. This is analogous to quarks, which 
are the elementary constituent particles of the baryons 
although they are strictly speaking always confined. 

Observing deconfined or almost deconfined spinons in 
experiments is in general difficult [6i]. In ID systems, 
e.g., the Heisenberg chain, it is well understood (based 
on the exact Bethe ansatz solution and numerical cal- 
culations 1^, III ) that spinons lead to a broad continuum 
in the dynamic spin structure factor S{q,uj). This con- 
tinuum has been observed in neutron scattering experi- 
ments on quasi-lD quantum antiferromagnets [8|. In 2D 
systems, there is no known reference model with decon- 
fined spinons in which S{q, w) can be computed exactly. 
One nevertheless expects a broad continuum also in this 
case, and such experimental signatures have been claimed 
in some quasi-2D systems 9]. The issue is complicated, 
however, by the fact that a continuum is also expected 
due to multi-magnon processes [lOJ. 

In this Letter, we discuss spinon detection in numeri- 
cal model calculations. This has also been a challenging 
problem, the solution of which will greatly help to un- 
derstand the conditions under which spinons can exist as 
independent elementary particles. Recently, signatures 



in thermodynamic properties were observed 11| in a 2D 
J-Q model (a spin- 1/2 Heisenberg model including four- 
spin interactions [l^]) at the point separating its Neel 
and VBS ground states. This model may, thus, exhibit 
the deconfined quantum-criticality proposed by Senthil 
et al. !5|. It is still desirable to have a more direct way 
to unambiguously (independently of any phenomenolog- 
ical ansatz or theory) detect spinons in numerical stud- 
ies of spin models (and eventually in doped systems). 
Here we introduce a method based on quantum Monte 
Carlo (QMC) simulations in the basis of valence bonds 
(singlet pairs) |l3| ■ generalized to include one or two un- 
paired spins |14l . |l5| . We show that an unpaired spin can 
constitute the core of a spinon wave packet, the size of 
which can be computed with our method. Analyzing the 
separation of two such wave packets we obtain quantita- 
tive information on the confinement or deconfinement of 
spinons in a magnetically disordered state. Importantly, 
our definitions also reproduce the expectation that the 
spinon should not be a low-energy particle in the ordered 
Neel state. 

Models — The primary model we use to test our method 
is the J-Qa chain, a ID member of the broad class of J-Q 
models introduced in Refs. [12|, [l6|. Defining a two-spin 
singlet projector dj = 1/4 — Si 



Sj, the hamiltonian is 
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where J, Qa > and we define g — Qs/ J. The ground 
state of this system is in the class of the standard critical 
Heisenberg chain for g < gc and is a doubly-degenerate 
VBS for g > gc- Using Lanczos diagonalization to extract 
the lowest singlet and triplet excitations and studying 
their crossings in the standard way for this kind of transi- 
tion (see, e.g., [17|), we obtain gc « 0.1645 (in agreement 



with a recent QMC study of the critical properties of the 
same model Il8|). We have also studied the J-Q2 model, 
i.e., using two singlet projectors in the Q-term in ([1]), for 
which gc ~ 0.84831. We here focus on the J-Q3 model 
because it is more strongly VBS ordered at J = 0. 




FIG. 1: (Color online) Spin (a) and dimer (b) correlations of 
systems with A'' = 1024 spins. Results for the J-Q3 model 
in the VBS phase ( J = 0, <? = 4, 1) and at criticality (^c) 
are shown along with the behavior in the Neel state of the 
long-range model with a = 3/2. The curves in (a) are fits to 
the form oc e"""'^ (with ^ ~ 4 at J = 0). The straight lines at 
the Qc data show the expected ~ 1/r critical behavior [27| . 

We also wish to study an ordered Neel state, which in 
an SU(2) invariant ID system can only be achieved with 
long-range interactions. The hamiltonian 
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Jr > 0, (2) 
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was studied in [l9|. With Jr — 1/r", a quantum phase 
transition from the critical state for a > Oc to a Neel 
state for a < ac was observed, with ac ~ 2.2. Here we 
use a slightly different model, with Jr = 1/r" for odd 
r but Jr — for even r, to make the system amenable 
to QMC simulations in the valence-bond basis [13]. We 
choose a = 3/2, for which the system is Neel ordered. 

To demonstrate the ground states of interest — VBS, 
critical, and Neel — in Fig. [T] we plot the spin and dimer 
correlation functions, defined by 

C(r) = (S, -S.+O, (3) 

i?(r)-((S, •S,+i)(S,+,, •S,+i+,)), (4) 

and computed using the QMC method discussed below. 
We multiply C{r) by (—1)'^ to cancel the signs of the 
correlations and graph (— l)'"[£'(r) — D{r + 1)], which for 
large r can be regarded as the VBS order parameter. 

QMC method — The valence-bond QMC algorithm and 
its generalizations to 5' > states have been discussed in 
several papers ISMlSi |2C| . Here we review key aspects of 
the basis and the form of the generated ground states. 
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FIG. 2: (Color online) Illustration of the basis for states with 
(a) 5* = (even A^), (b) S = 1/2 (odd iV), and (c) 5 = 1 (even 
N) . The bonds and unpaired spins of the bra and ket states 
are shown below and above the line of sites, respectively. 

Acting with a high power of the hamiltonian H™' on a 
trial state |^t), with H written as a sum of singlet projec- 
tors (individual ones and products of three, for J and Q 
interactions, respectively), the ground-state normaliza- 
tion (^ol^o) is sampled (for m large enough for H"'-\'$t) 
to be completely dominated by |^o))- In an S* = state 
for even N, the states are expressed as superpositions of 
bipartite valence-bond states \Va), i.e., products of N/2 
singlets (a, b) = {taib — ibta)/V2 where a and b are sites 
on sublattice A and B, respectively. We use trial states 
of the amplitude-product form pii]. 

The valence-bond basis is non-orthogonal, and the nor- 
malization of the projected ground state is therefore of 
the form (*o|*o) = Y.a0 fpfaiVplVa), where fpja are 
not known explicitly. Implicitly, the probability of gen- 
erating a pair of states is P(T4,, V^) = fpfaiVplVa). The 
overlap (VfilVa) = 2^°'^/"^, where iVo is the number of 
loops in the transition graph of the two states. Fig. HJa) 
shows a case with A^o = 1 • Matrix elements of the form 
(1^|A|T4.) for many observables A of interest depend on 
the loop structure of the transition graph 2lJ, |22 |. 

For S > and magnetization rriz = S the states have 
2mz unpaired f spins and {N — 2mz)/2 singlet bonds (as 



discussed, e.g., in [ij, ll5|). For odd N, which we use 



for S — 1/2, the system is in principle frustrated by pe- 
riodic boundaries. This is a finite-size effect, however, 
which vanishes when N -^ 00 {a.t least for observables 
probing distances r ^ N). The QMC loop updates [20| 
automatically exclude frustrated negative-sign configu- 
rations, and this should, thus, be the most rapid way to 
approach N — 00. Configurations for S* = 1/2 and 5 = 1 
states are illustrated in Fig. [2][b,c). We note that the 
valence bond basis with two unpaired spins was used in 
a pioneering variational study on spinon deconfinement 
in a VBS state of a ID frustrated model [l]. 

Spinon statistics — The first aspect of our method re- 
lies on the representation of 5* = 1/2 states in terms of 
valence-bond states with an unpaired spin [15]. One can 
determine whether there is a well-defined wave packet 
(localizable particle) carrying the spin. The second as- 
pect is to characterize the correlations of two spinous in 
an 5 = 1 state, to determine whether they are confined. 
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FIG. 3: (Color online) Overlap P{r) = {V'i/2(i+r-)|Vi/2(i)) for 

(a) different VBS states of the J-Q3 model of size A*' — 1025, 

(b) at gc for different A'^, and (c) in the Neel state of the long- 
range model (a = 3/2) for different A'". The curves in (a) are 
fits to oc e~^' (with A « 9 at J = 0) and the line in (b) 
shows the form oc l/\/r- 

and, if so, to extract the size of the bound state. 

The S = 1/2 ground state (with momentum k = 0) 
can be written as 1*1/2) = XIJV'1/2 (?")), where r is the 
location of the unpaired spin [15| . Denoting a basis state 
with the spinon at r as |VQ.(r)), we have |'0i/2('')) = 
J2a fr\^a{r)) and the overlap of two states with different 
location of their spinon cores is 

(V'i/2(/)l^i/2M) ^Y.fr'fr{Vf3ir')\V^ir)). (5) 

What we propose is that this quantity allows for a generic 
way to test whether a spinon is a well defined particle. 
Such a particle should have a finite wave packet (i.e., a 
minimum size of a region to which the S = 1/2 degree 
of freedom can be confined) , which typically should lead 
to an exponential decay of the overlap with the separa- 
tion \r' — r\ (with a power-law decay corresponding to a 
marginal case). This follows in a VBS state because the 
basis-state overlap {Vfs{r')\Va{r)) is dictated by the num- 
ber of loops in the transition graph. An S = 1/2 transi- 
tion graph has a string of bonds terminating in unpaired 
spins [15|, as seen in^b). In a VBS state, the loops are 
typically short, and the presence of a string will reduce 



the number of loops in proportion to the length of the 
string, and, thus, {Vf3{r')\Va{r)) and (O, should decay 
exponentially with the separation \r' — r\. One can then 
also expect a power-law decay in a critical VBS state. 

The overlap ([5]) can be computed by accumulating the 
distribution P{r) of separations r of the unpaired spins in 
the S ~ 1/2 transition graphs. The above expected be- 
haviors are indeed realized in the J-Q3 model, as shown 
in Fig. [21 In VBS states for large N, the overlap vanishes 
for odd distances, implying that the bra and ket spinous 
are on the same sublattice in the infinite system. For the 
even distances the overlap is of the form P{r) ex e^"^/^, 
and A is essentially the size of an exponentially decay- 
ing wave packet. The size A is roughly twice the spin 
correlation length in the cases we have studied. 

In the critical state, the overlap has the form P{r) ^ 
l/\/r and the wave packet is only marginally defined. 
The Heisenberg chain is known to have spinon excitations 
[2] and, thus, it appears that one can still consider such 
a broad algebraic wave packet as a particle. The total 
weight of all odd-r overlaps is roughly constant, « 1/4. 

In the Neel state P{r) is almost flat and even and odd 
r have almost the same weight. The unpaired spin in the 
Neel state is, thus, not localizable within a wave packet, 
in agreement with the expectation that the spinon should 
not be an elementary excitation of this state. The un- 
paired spin is strongly aligned with the Neel order of the 
rest of the system [23j and cannot be regarded as an in- 
dependent spatial S = 1/2 degree of freedom. 

For an S* = 1 state with two unpaired spins we have 

{Mr'A^r'i^)\MrA,rB)) 
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Jr.rs{M^'AyB)\Vo.{rA,rB)). (6) 



where, as indicated, in both the bra and the ket state 
one spinon is on sublattice A and one on B. Here we 
can define several probability distributions depending on 
a single distance, e.g., |r^ — rsj or |r^ — f'g\, integrat- 
ing over the remaining two free-spin locations. To inves- 
tigate the confinement length we define PAB(r) as the 
average of the distributions of the above two distances. 
The results, shown in Fig. 21 indicate deconfined spinous 
with weak mutual repulsion, which makes the distribu- 
tion broadly peaked at r = N/2. Size-scaled distributions 
NPab {r) for different N fall almost on top of each other 
when graphed versus r/N. For confined spinous, the con- 
finement length will be refiected in an asymptotic decay 
PAB{r) ^ e~''/'^, where A is the size of the bound state. 
The two length-scales we have discussed — the size of 
the spinon wave packet A and the bound state A (for 
confined spinous) — are also visible in the z-component 
spin correlation function Cz{r) = (SfSf^j.). We demon- 
strate this for both S = 1/2 and S = 1 states in Fig. |5l 
The reason for the spinon contributions can be under- 
stood from Figs.JH^bjC). The unpaired spins in the S — 1 
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FIG. 4: (Color online) Size-scaled distribution of the distance 
between two spinons in the S = 1 VBS state at (; = 4. 

state (c) dominate the long-distance correlation function 
if the confinement length is larger than the correlation 
length of the backgromid VBS. In addition, the short- 
distance correlations of both the 5 = 1/2 and 5 = 1 
states are modified by the presence of strings. Indeed, 
as we demonstrate in Fig. [SI by subtracting off Cz {r) of 
the 5 = state, the remaining short-distance correlations 
contain an exponentially decaying contribution which for 
both 5 = 1/2 and 5 = 1 is roughly twice the correlation 
length, i.e., similar to the wave packet size A. In the 
5=1 state, the correlation function remains non-zero, 
ex 1/iV, as r — )■ oo, reflecting deconfined spinons. Note 
that there is a change in phase of Cz(r), at some r which 
is related to N and A (and hence depends on g). 

Conclusions and discussion — We have introduced a 
method to determine whether a spinon is a well-defined 
emergent particle (excitation) of a quantum spin system, 
and, if so, whether two spinons in an 5 = 1 excitation 
are deconfined or form a bound state (the size of which 
can be computed). The discussion was framed around the 
valence-bond basis and QMC simulations with it, but the 
definitions are independent of this basis. Our arguments 
only rely on the fact that one can write a state for, e.g., 
5 = 1/2 as ^^ \i^o{f))®\ tr) (for momentum fc = 0, with 
self-evident generalization to fc 7^ 0), where \ipQ{r)) is an 
5 = state of all spins except the one at r (and simi- 
lar decompositions for higher 5). One can, thus, com- 
pute the quantities we have investigated here with other 
methods as well. The crucial observation is that states 
IV'o(^)) ® I tr) for different r are non-orthogonal. If the 
unpaired spin fr is localized within a spinon wave packet 
(by definition for a spinon) , then the overlaps give direct 
information on the size of this wave packet. The spinon 
is not an independent particle if the wave packet is uni- 
formly delocalized over the whole system as A^ -^ 00, as 
we have demonstrated here for a Neel state. 

Our method does not rely on any knowledge or theory 
of the nature of the spinon (other than it carrying spin 
5 = 1/2). The wave-function overlaps ([5]) and dU are 
completely general and applicable to any system in any 
number of dimensions. For ID systems there are alterna- 
tive methods to study spinons using the fact that they are 
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FIG. 5: (Color online) Spin correlations in an A'' = 513 system 
with a single spinon {S = 1/2) and in A'' = 512 systems with 
two spinons (S = 1) in the J = VBS state. The N = 512, 
5 = correlation has been subtracted off to isolate the spinon 
contributions. For S = 1 there is a phase change at r « 42. 

kink and antikink solitons [2J], which can be created by 
boundary conditions. The spinon wave function, which 
is similar to that of a particle in a box [25], does not, 
however, contain any direct information on the intrinsic 
size of the spinon "particle" . A criterion of deconfine- 
ment based on impurity (un)binding was also presented 
recently [26|, but that approach cannot unambiguously 
determine whether a spinon is a well defined particle. 
Our approach also avoids potential differences between 
spinon-spinon and spinon-impurity affinities 
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